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Abstract: Results are presented of a search for new particles decaying to large numbers 
of jets in association with missing transverse momentum, using 4.7 fb -1 of pp collision 
data at y/s = 7 TeV collected by the ATLAS experiment at the Large Hadron Collider in 
2011. The event selection requires missing transverse momentum, no isolated electrons or 
muons, and from >6 to >9 jets. No evidence is found for physics beyond the Standard 
Model. The results are interpreted in the context of a MSUGRA/CMSSM supersymmetric 
model, where, for large universal scalar mass mo, gluino masses smaller than 840 GeV 
are excluded at the 95% confidence level, extending previously published limits. Within 
a simplified model containing only a gluino octet and a neutralino, gluino masses smaller 
than 870 GeV are similarly excluded for neutralino masses below 100 GeV. 
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1 Introduction 

Many extensions of the Standard Model of particle physics predict the presence of TeV- 
scale strongly interacting particles that decay to lighter, weakly interacting descendants. 
Any such weakly interacting particles that are massive and stable can contribute to the 
dark matter content of the universe. The strongly interacting parents would be produced in 
the proton-proton interactions at the Large Hadron Collider (LHC), and such events would 
be characterized by significant missing transverse momentum E!f ms from the unobserved 
weakly interacting daughters, and jets from emissions of quarks and/or gluons. 

In the context of i?-parity conserving [1-5] supersymmetry [5-10], the strongly inter- 
acting parent particles are the squarks q and gluinos g, they are produced in pairs, and the 
lightest supersymmetric particles can provide the stable dark matter candidates [11, 12]. 
Jets are produced from a variety of sources: from quark emission in supersymmetric cas- 
cade decays, production of heavy Standard Model particles (W, Z or t) which then decay 
hadronically, or from QCD radiation. Examples of particular phenomenological interest 
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include models where squarks are significantly heavier than gluinos. In such models the 
gluino pair production and decay process 



can dominate, producing large jet multiplicities when the resulting top quarks decay 
hadronically. In the context of MSUGRA/CMSSM models, a variety of different cascade 
decays, including the gg initiated process above, can lead to large jet multiplicities. 

A previous ATLAS search in high jet multiplicity final states [13] examined data taken 
during the first half of 2011, corresponding to an integrated luminosity of 1.34 fb. This 
paper extends the analysis to the complete ATLAS 2011 pp data set, corresponding to 
4.7 fb -1 , and includes improvements in the analysis and event selection that further increase 
sensitivity to models of interest. 

Events are selected with large jet multiplicities ranging from > 6 to > 9 jets, in associ- 
ation with significant i?™ lss . Events containing high transverse momentum (pr) electrons 
or muons are vetoed in order to reduce backgrounds from (semi-leptonically) decaying top 
quarks or W bosons. Other complementary searches have been performed by the ATLAS 
collaboration in final states with E!f lss and one or more leptons [14, 15]. Further searches 
have been carried out by ATLAS using events with at least two, three or four jets [16], 
or with at least two 6-tagged jets [17]. Searches have also been performed by the CMS 
collaboration, including a recent analysis in fully hadronic final states [18]. 

2 The ATLAS detector and data samples 

The ATLAS experiment [19] is a multi-purpose particle physics detector with a forward- 
backward symmetric cylindrical geometry and nearly Air coverage in solid angle. 1 The lay- 
out of the detector is dominated by four superconducting magnet systems, which comprise 
a thin solenoid surrounding inner tracking detectors, and a barrel and two end-cap toroids 
supporting a large muon spectrometer. The calorimeters are of particular importance to 
this analysis. In the pseudorapidity region \rj\ < 3.2, high-granularity liquid-argon (LAr) 
electromagnetic (EM) sampling calorimeters are used. An iron/scintillator-tile calorime- 
ter provides hadronic coverage for |r/| < 1.7. The end-cap and forward regions, spanning 
1.5 < \t]\ < 4.9, are instrumented with LAr calorimetry for both EM and hadronic mea- 
surements. 

The data sample used in this analysis was taken during April - October 2011 with 
the LHC operating at a proton-proton centre-of-mass energy of -y/s = 7 TeV. Application 
of beam, detector and data-quality requirements resulted in a corresponding integrated 
luminosity of 4.7±0.2fb _1 [20]. The analysis makes use of dedicated multi-jet triggers that 
required either at least four jets with pi > 45 GeV or at least five jets with > 30 GeV, 

1 ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point in the 
centre of the detector and the z-axis along the beam pipe. Cylindrical coordinates (r, <j>) are used in the 
transverse plane, <f) being the azimuthal angle around the beam pipe. The pseudorapidity r\ is defined in 
terms of the polar angle 9 by r\ — — lntan(#/2). 
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where the energy is measured at the electromagnetic scale 2 and the jets must have \r]\ < 3.2. 
In all cases the trigger efficiency was greater than 98% for events satisfying the offline jet 
multiplicity selections described in Section 4. 

3 Object reconstruction 

The jet, lepton and missing transverse momentum definitions are based closely on those of 
Ref. [13], with small updates to account for evolving accelerator and detector conditions. 

Jet candidates are reconstructed using the anti-kt jet clustering algorithm [21, 22] 
with radius parameter of 0.4. The inputs to this algorithm are clusters of calorimeter 
cells seeded by cells with energy significantly above the noise level. Jet momenta are 
reconstructed by performing a four- vector sum over these topological clusters of calorimeter 
cells, treating each as an (E,p) four- vector with zero mass. The jet energies are corrected 
for the effects of calorimeter non-compensation and inhomogeneities by using pp- an d 
77-dependent calibration factors based on Monte Carlo (MC) simulations validated with 
extensive test-beam and collision-data studies [23]. Only jet candidates with pp > 20 GeV 
and \rj\ < 4.9 are retained. Further corrections are applied to any jet falling in problematic 
areas of the calorimeter. The event is rejected if, for any jet, this additional correction 
leads to a contribution to E™ 1SS that is greater than both 10 GeV and 0.1 E™ 1SS . These 
criteria, along with selections against non-collision background and calorimeter noise, lead 
to a loss of signal efficiency of ~8% for the models considered. When identification of 
jets containing heavy flavour quarks is required, either to make measurements in control 
regions or for cross checks, a tagging algorithm exploiting both impact parameter and 
secondary vertex information is used. Jets are tagged for \r/\ < 2.5 and the parameters of 
the algorithm are chosen such that 70% of 6-jets and ~1% of light flavour or gluon jets, 
are selected in ti events in Monte Carlo simulation [24]. Jets initiated by charm quarks are 
tagged with about 20% efficiency. 

Electron candidates are required to have pp > 20 GeV and \rj\ < 2.47, and to sat- 
isfy the 'medium' electron shower shape and track selection criteria of Ref. [14]. Muon 
candidates are required to have pp > 10 GeV and |?7| < 2.4. Additional requirements are 
applied to muons when defining leptonic control regions. In this case muons must have 
longitudinal and transverse impact parameters within 1mm and 0.2 mm of the primary 
vertex, respectively, and the sum of the transverse momenta of other tracks within a cone 
of AR = 0.2 around the muon must be less than 1.8 GeV, where AR = \J (Arj) 2 + (A<p) 2 . 

The measurement of the missing transverse momentum two- vector p^? lss and its magni- 
tude (conventionally denoted E™ 1SS ) is then based on the transverse momenta of all electron 
and muon candidates, all jets with \rj\ < 4.5 which are not also electron candidates, and 
all calorimeter clusters with \tj\ < 4.5 not associated to such objects [25]. 

2 The electromagnetic scale is the basic calorimeter signal scale for the ATLAS calorimeters. It has 
been established using test-beam measurements for electrons and muons to give the correct response for 
the energy deposited in electromagnetic showers, although it does not correct for the lower response of the 
calorimeter to hadrons. 
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Signal region 


7j55 


8j55 


9j55 


6j80 


7j80 


8j80 


Number of isolated leptons (e, //) 


= 


Jet pt 


> 55 GeV 


> 80 GeV 


Jet 77 


<2.8 


Number of jets 


>7 


>8 


>9 


>6 


>7 


>8 




> 4 GeV 1 / 2 



Table 1. Definitions of the six signal regions. 



Following the steps above, overlaps between candidate jets with |r/| < 2.8 and leptons 
are resolved as follows. First, any such jet candidate lying within a distance AR = 0.2 of an 
electron is discarded, then any lepton candidate remaining within a distance AR = 0.4 of 
such a jet candidate is discarded. Thereafter, all jet candidates with \r]\ > 2.8 are discarded, 
and the remaining electron, muon and jet candidates are retained as reconstructed objects. 

4 Event selection 

Following the object reconstruction described in Section 3, events are discarded if they 
contain any jet failing quality criteria designed to suppress detector noise and non-collision 
backgrounds, or if they lack a reconstructed primary vertex with five or more associated 
tracks. 

For events containing no isolated electrons or muons, six non-exclusive signal regions 
(SRs) are defined as shown in Table 1. The first three require at least seven, eight or nine 
jets, respectively, with > 55 GeV; the latter three require at least six, seven or eight 
jets, respectively, with j»t > 80 GeV. The final selection variable is E<f lss /\/Hic, the ratio 
of the magnitude of the missing transverse momentum to the square root of the scalar 
sum Ht of the transverse momenta of all jets with > 40 GeV and \r]\ < 2.8. This 
ratio is closely related to the significance of the missing transverse momentum relative to 
the resolution due to stochastic variations in the measured jet energies [25]. The value of 
E™ 1ss /\/Ht is required to be larger than 4 GeV 1//2 for all signal regions. 

A previous ATLAS analysis of similar final states [13] required jets to be separated 
by AR > 0.6 to ensure that the trigger efficiency was on its plateau. It has since been 
demonstrated that the requirement of an offline jet multiplicity at least one larger than 
that used in the trigger is sufficient to achieve a 98% trigger efficiency. Investigations on 
the enlarged data sample, in comparison to the previous incarnation of the strategy used 
here, allow various improvements to be made; in particular, the requirement on jet-jet 
separation is modified so as to increase the acceptance for signal models of interest by a 
factor two to five, without introducing any significant trigger inefficiency. 

The dominant backgrounds are multi-jet production, including purely strong interac- 
tion processes and fully hadronic decays of tt; semi- and fully-leptonic decays of ti; and 
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leptonically decaying W or Z bosons produced in association with jets. Non-fully-hadronic 
ti, and W and Z are collectively referred to as 'leptonic' backgrounds. Contributions from 
gauge boson pair and single top quark production are negligible. The determination of the 
multi-jet and 'leptonic' backgrounds is described in Sections 6 and 7, respectively. 

5 Monte Carlo simulations 

Monte Carlo simulations are used as part of the 'leptonic' background determination pro- 
cess, and to assess sensitivity to specific SUSY signal models. The 'leptonic' backgrounds 
are generated using Alpgen2.13 [26] with the PDF set CTEQ6L1 [27]. Fully-leptonic ti 
events are generated with up to five additional partons in the matrix element, while semi- 
leptonic ti events are generated with up to three additional partons in the matrix element. 
W + jets and Z — > vv + jets are generated with up to six additional partons, and the 
Z — > l + £~ + jets (for £ £ {e, //, r}) process is generated with up to five additional partons 
in the matrix element. In all cases, additional jets are generated via parton showering, 
which, together with fragmentation and hadronization, is performed by Herwig [28, 29]. 
Jimmy [30] is used to simulate the underlying event. The W + jets, Z + jets and ti 
backgrounds are normalized according to their inclusive theoretical cross sections [31, 32]. 
The estimation of the 'leptonic' backgrounds in the signal regions is described in detail in 
Section 7. 

Supersymmetric production processes are generated using Herwig++2.4.2 [33]. Sig- 
nal cross sections are calculated to next-to-leading order in the strong coupling constant 
as, including the resummation of soft gluon emission at next-to- leading-logarithmic ac- 
curacy (NLO+NLL) [34-38]. 3 An envelope of cross-section predictions is defined using 
the 68% confidence-level (CL) ranges of the CTEQ6 . 6 [39] (including the as uncertainty) 
and MSTW2008 NLO [40] PDF sets, together with independent variations of the factor- 
ization and renormalization scales by factors of two or one half. The nominal cross- 
section value is then taken to be the midpoint of the envelope, and the uncertainty as- 
signed is half the full width of the envelope, following closely the PDF4LHC recommen- 
dations [41]. MSUGRA/CMSSM particle spectra and decay modes are calculated with 
ISAJET++7.75 [42]. For illustrative purposes, plots of kinematic quantities show the distri- 
bution expected for an example MSUGRA/CMSSM point that has not been excluded in 
previous searches. This reference point is defined by 4 : mo = 2960 GeV, m 1 / 2 = 240 GeV, 
A = 0, tan/3 = 10, and fi > 0. 

3 The NLL correction is used for squark and gluino production when the average of the squark masses 
in the first two generations and the gluino mass lie between 200 GeV and 2 TeV. In the case of gluino-pair 
(associated squark-gluino) production processes, the calculations were extended up to squark masses of 4.5 
TeV (3.5 TeV). For masses outside this range and for other types of production processes (i.e. electroweak 
and associated strong and electroweak), cross sections at NLO accuracy obtained with Prospino2.1 [34] 
are used. 

4 A particular MSUGRA/CMSSM model point is specified by five parameters: the universal scalar mass 
mo, the universal gaugino mass m^ 2 , the universal trilinear scalar coupling Ao, the ratio of the vacuum 
expectation values of the two Higgs fields tan/?, and the sign of the higgsino mass parameter jj,. 
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Figure 1. E<f lss / \J distributions in example multi-jet control regions, (a) For exactly six 
jets with pt > 55 GeV, compared to a prediction based on the E™ lss / \J Ht distribution for exactly 
five jets with > 55 GeV. (b) For exactly five jets with px > 80 GeV, compared to a prediction 
based on four jets with pt > 80 GeV. The multi-jet predictions have been normalized to the data 
in the region E^ lss /y/Hr < 1.5 GeV 1 / 2 after subtraction of the predicted 'leptonic' backgrounds. 
The most important 'leptonic' backgrounds are also shown, based on MC simulations. Variable bin 
sizes are used with bin widths (in units of GeV 1//2 ) of 0.25 (up to 4), 0.5 (from 4 to 5), 1 (from 5 
to 6), and then 2 thereafter. The error bars on the data points show the Poisson coverage interval 
corresponding to the number of data events observed in each bin. 



All Monte Carlo samples employ a detector simulation [43] based on GEANT4 [44] and 
are reconstructed with the same algorithms as the data. 

6 Multi-jet backgrounds 

The dominant background at intermediate values of E!f lss is multi-jet production including 
purely strong interaction processes and fully hadronic decays of ti. These processes are not 
reliably predicted with existing Monte Carlo calculations, and so their contributions must 
be determined from collision data. Indeed, the selection cuts have been designed such that 
multi-jet processes can be determined reliably from supporting measurements. 

The method for determining the multi-jet background from data is motivated by the 
following considerations. In events dominated by jet activity, including hadronic decays of 
top quarks and gauge bosons, the E™ ss resolution is approximately proportional to \/Ht, 
and is almost independent of the jet multiplicity. The distribution of the ratio E^ 1SS /^/I^^- 
has a shape that is almost invariant under changes in the jet multiplicity, as shown in 
Figure 1. The multi-jet backgrounds therefore can be determined using control regions 
with lower E™ 1ss /\/Ht and/or lower jet multiplicity than the signal regions. 5 The control 

5 Residual variations in the shape of the E™' BB / \/Ht are later used to quantify the systematic uncertainty 
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regions are assumed to be dominated by Standard Model processes, an assumption that 
is corroborated by the agreement of multi-jet cross section measurements with up to six 
jets [45] with Standard Model predictions. 

As an example, the estimation of the background expected in the 8j55 signal region 
is obtained as follows. A template describing the shape of the i?™ lss /\/i?T distribution is 
obtained from those events that contain exactly six jets, using the same 55 GeV pp threshold 
as the target signal region. That six-jet E™ 1SS /yHr template is normalized to the number 
of eight-jet events observed in the region 

£2p is 7V5r < 1.5 GeV 1/2 after subtraction of the 
'leptonic' background expectation. The normalized template then provides a prediction for 
the multi-jet background for the 8j55 signal region for which E™ 1s */^Ht > 4 GeV 1 / 2 . 

A similar procedure is used for each of the signal regions, and can be summarized as 
follows. For each jet pr threshold p< G {55 GeV, 80 GeV}, control regions are defined for 
different numbers n ]Ct of jets found above p<. The number of events A^ <)njet (s TO i n , s max ) 
for which E™ ss / ^fEE? (in units of GeV 1 / 2 ) lies between s m i n and s max is determined, and 
the predicted 'leptonic' contributions L p<>n . et (s m ; n , s max ) subtracted 

-^j^onjet (^min) ^max) — -^p<,nj e t (^mini ^max) ^- / P<,«j c t (^min> S max ). 
Transfer factors 

JV&, Wjct (4, oo) 
N p ^ n}ct (0, 1.5) 

connect regions with the same p < and nj e t with different E™ 1SS / \/St- The multi-jet pre- 
diction for the signal region is found from the product of the T p< iTl . et , with the same p < as 
the signal region and nj e t = 6 when p < = 55 GeV (rij et = 5 when p < = 80 GeV) times the 
number of events (after subtracting the expected contribution from 'leptonic' background 
sources) satisfying signal region jet multiplicity requirements but with £™ ss /\/5r < 
1.5 GeV 1 / 2 . 

6.1 Systematic uncertainties on multi-jet backgrounds 

The method is validated by determining the accuracy of predictions for regions with jet 
multiplicities and/or E™ ss /\/H^ smaller than those chosen for the signal regions. Figure 1 
shows that the shape of the E™ ss /\/St distribution for p< = 55 GeV and nj e t = 6 is 
predicted to an accuracy of better than 20% from that measured using a template with the 
same value of p < and n- ]e t = 5. Similarly, the distribution for p< = 80 GeV and nj et = 5 
can be predicted for all E!f lss / V-f^r using a template with n- ]e t = 4. The templates are 
normalized for E™ 1SS / \/H^ < 1.5 GeV 1 / 2 , and continue to provide a good prediction of the 
distribution out to values of E^ 1SS /^/^^^' of 4 GeV 1 / 2 and beyond. Additional validation 
regions are defined for each p < and for jet multiplicity requirements equal to those of 
the signal regions, but for the intermediate values of (s m i n , s max ) of (1.5, 2), (2, 2.5) and 
(2.5, 3.5). Residual inaccuracies in the predictions are used to quantify the systematic 
uncertainty from the closure of the method. Those uncertainties are in the range 15%- 
25%, depending on p< and E™ ss /\/H^. 

associated with the method, as described in Section 6.1. 
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The mean number of proton-proton interactions per bunch crossing (//) increased dur- 
ing the 2011 run, reaching (fj) = 16 at the start of proton fills for runs late in the year. 
Sensitivity to those additional interactions is studied by considering the jet multiplicity as 
a function of (//), and of the number of reconstructed primary vertices. The consistency of 
the high-px tracks within the selected jets with a common primary vertex is also investi- 
gated. The effect of additional jets from pile-up interactions is found to be significant for 
1ow-pt jets but small for jets with pr > 45 GeV, and negligible for the jet selection used 
for the SRs. 

The presence of multiple in-time and out-of-time pp interactions also leads to a small 
but significant deterioration of the E™ 1SS resolution. The effectiveness of the E™ ss / \J 
template method described above is tested separately for subsets of the data with differ- 
ent values of the instantaneous luminosity, and hence of (fi). Good agreement is found 
separately for each subset of the data. Since the data set used to form the template has 
the same pile- up conditions as that used to form the signal regions, the changing shape of 
the E™ 1SS resolution is included in the data-driven determination and does not lead to any 
additional systematic uncertainty. 

Due to the presence of neutrinos produced in the decay of hadrons containing bottom 
or charm quarks, events with heavy-flavour jets exhibit a different E™ ms distribution. To 
quantify the systematic uncertainty associated with this difference, separate templates are 
defined for events with at least one 6-tagged jet and for those with none. The sum of 
the predictions for events with and without 6-tagged jets is compared to the flavour-blind 
approach, and the difference is used to characterize the systematic uncertainty from heavy 
flavour (10%-20%). Other systematic uncertainties account for imperfect knowledge of: 
the subtracted 'leptonic' contributions (10%), the potential trigger inefficiency (2%), and 
imperfect response of the calorimeter in problematic areas (1%). 

The backgrounds from multi-jet processes are cross checked using another data-driven 
technique [16] which smears the energies of individual jets from low-E™ ss multi-jet 'seed' 
events in data. Separate smearing functions are defined for 6-tagged and non-6-tagged jets, 
with each modelling both the Gaussian core and the non-Gaussian tail of the jet response, 
including the loss of energy from unobserved neutrinos. The jet smearing functions are 
derived from GEANT4 [44] simulations [43] . The Gaussian core of the function is tuned to 
di-jet data, and the non-Gaussian tails are verified with data in three-jet control regions 
in which the p™ lss can be associated with the fluctuation of a particular jet. There is 
agreement within uncertainties between the background predicted by this jet-smearing 
method and the primary method based on the shape invariance of E!f lss /\/H^. 

7 'Leptonic' backgrounds 

Non-fully-hadronic (i.e. semi-leptonic or di-leptonic) tt, and W and Z production are col- 
lectively referred to as 'leptonic' backgrounds. The process Z — > vv + jets contributes to 
the signal regions since it produces jets in association with E™ 1SS . Leptonic tt and W decays 
contribute to the signal regions when hadronic r decays allow them to evade the lepton 
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Table 2. Definitions of the validation regions and control regions for the 'leptonic' backgrounds: 
tt + jets, W + jets and Z + jets. The validation regions VR are defined by the first five selection 
requirements. A long dash ' — ' indicates that no requirement is made. The control regions CR 
differ from the VR in their treatment of the muons, and by having additional requirements on jets 
and E^ 1ss /^/Ht, as shown in the final two rows. 



veto, with smaller contributions from events in which electrons or muons are produced but 
are not reconstructed. 

The 'leptonic' background predictions employ the Monte Carlo simulations described 
in Section 5. To reduce uncertainties from Monte Carlo modelling and detector response, 
it is desirable to normalize the background predictions to data using control regions (CR) 
and cross-check them against data in other validation regions (VR). These control regions 
and validation regions are designed to be distinct from, but kinematically close to, the 
signal regions. Each is designed to provide enhanced sensitivity to a particular background 
process. 

The control and validation regions are defined as shown in Table 2. By using control 
regions that are kinematically similar to the signal regions, theoretical uncertainties, includ- 
ing those arising from the use of a leading-order (LO) generator, are reduced. The tt + jets 
and W + jets validation regions each require a single muon and no electrons. For the tt 
process the single-muon selection is primarily sensitive to the semi-leptonic decay. 6 The 
tt + jets validation region is further enhanced by the requirement of at least one 6-tagged 
jet, whereas for W + jets enhancement a 6-tag veto is applied. Since it is dominantly 
through hadronic r decays that W and tt contribute to the signal regions, the correspond- 
ing control regions are created by recasting the muon as a (r-)jet. For Z — )■ vv + jets 
the validation regions select events from the closely related process Z — > fifi + jets. The 
related control regions are formed from these validation regions by recasting the muons as 
neutrinos. 

In detail, for those control regions where the Monte Carlo simulations predict at least 

6 The procedure is also sensitive to those di-leptonic ti decays in which one lepton was not observed 
in the VR. After the VR — > CR replacement (fi — > jet), the procedure captures the leading di-leptonic ti 
contributions to the SR. 
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one event for 4.7 fb -1 , the leptonic background prediction s, for each signal region from 
each background is calculated by multiplying the number of data events c^ ata found in the 
corresponding control region by a Monte Carlo-based factor $ 



,MC 

i 



. _ data .MC 

"t — "% 

This transfer factor is defined to be the ratio of the number of MC events found in the 
signal region to the number of MC events found in the control region 

„MC 
^MC _ s i 



°i 

In each case, the event counts are corrected for the expected contamination by the other 
background processes. Whenever less than one event is predicted in the control region, 
the Monte Carlo prediction for the corresponding signal region is used directly, without 
invoking a transfer factor. 

For the ti + jets background, the validation region requires exactly one isolated muon, 
at least one 6-tagged jet, and no selected electrons. The transverse mass for the muon trans- 
verse momentum p£ and the missing transverse momentum two- vector p™ lss is calculated 
using massless two- vectors 

m 2 = 2\p»\\p^ iss \ - 2p% ■ pf iss , 

and must satisfy 50 GeV < wit < 100 GeV. Figure 2 shows the jet multiplicity in the ti 
validation regions, and it is demonstrated that the Monte Carlo provides a good description 
of the data. 

The ti control regions used to calculate the background expectation differ from the 
validation regions as follows. Since the dominant source of background is from hadronic 
r decays in the control regions, the muon is used to mimic a jet, as follows. If the muon 
has sufficient pr to pass the jet selection threshold p<, the jet multiplicity is incremented 
by one. If the muon p^ is larger than 40 GeV it is added to flr- The selection variable 
E™ ss /i/5t is then recalculated, and required to be larger than the threshold value of 
4 GeV 1//2 . Distributions of the jet multiplicity in the ti control regions may also be found 
in Figure 2. 

The W + jets validation regions and control regions are defined in a similar manner 
to those for ti + jets, except that a 6-jet veto is used rather than a 6-jet requirement (see 
Table 2). Figure 3 shows that the resulting jet multiplicity distributions are well described 
by the Monte Carlo simulations. 

The Z + jets validation regions are defined (as shown in Table 2) requiring precisely 
two muons with invariant mass consistent with mz- The dominant backgrounds from 
Z + jets arise from decays to neutrinos, so in forming the Z + jets control regions from the 
validation regions, the vector sum of the pr of the muons is added to the measured p™ ss , 
to model the E!f lss expected from Z — > vv events. The selection variable E^ 1SS /^ I /H^ 
is then recalculated and required to be greater than 4 GeV 1 / 2 for events in the control 
region. Figure 4 shows that the resulting jet multiplicity distributions in both validation 
and control regions are well described by the Monte Carlo simulations. 
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For each of the 'leptonic' backgrounds further comparisons are made between Monte 
Carlo and data using the lower jet threshold of 45 GeV, showing agreement within 
uncertainties for all multiplicities (up to nine jets for tt, see Figures 2(a) and 2(b). The 
Alpgen Monte Carlo predictions for Z + jets and W + jets were determined with six 
additional partons in the matrix element calculation, and cross checked with a calculation in 
which only five additional partons were produced in the matrix element - in each case with 
additional jets being produced in the parton shower. The two predictions are consistent 
with each other and with the data, providing further supporting evidence that the parton 
shower offers a sufficiently accurate description of the additional jets. 

7.1 Systematic uncertainties on 'leptonic' backgrounds 

The use of control regions is effective in reducing uncertainties from Monte Carlo modelling 
and detector response. When predictions are taken directly from the Monte Carlo, the 
'leptonic' background determinations are subject to systematic uncertainties from Monte 
Carlo modelling of: the jet energy scale (JES, 40%), the jet energy resolution (JER, 4%), 
the number of multiple proton-proton interactions (3%), the 6-tagging efficiency (5% for 
tt), the muon trigger and reconstruction efficiency and the muon momentum scale. The 
numbers in parentheses indicate the typical values of the SR event yield uncertainties prior 
to the partial cancellations that result from the use of control regions. 

The JES and JER uncertainties are calculated using a combination of data-driven and 
Monte Carlo techniques [23], using the complete 2011 ATLAS data set. The calculation 
accounts for the variation in the uncertainty with jet and rj, and that due to nearby 
jets. The Monte Carlo simulations model the multiple proton-proton interactions with a 
varying value of (fj) which is well matched to that in the data. The residual uncertainty 
from pile-up interactions is determined by reweighting the Monte Carlo samples so that (/i) 
is increased or decreased by 10%. The uncertainty in the integrated luminosity is 3.9% [20]. 

When transfer factors are used to connect control regions to signal regions, the effects 
of these uncertainties largely cancel in the ratio. For example, the impact of the jet energy 
scale uncertainty is reduced to « 6%. 

8 Results, interpretation and limits 

Figure 5 shows the E<f 1SS /^/I^^ distributions after applying the jet selections for the six 
different signal regions (see Table 1) prior to the final Ef iBB /^/H^ > 4 GeV 1/2 requirement. 
Figure 6 shows the jet multiplicity distributions for the two different jet pt thresholds 
after the final i?™ ss /\/ffr requirement. It should be noted that the signal regions are not 
exclusive. For example, in Figure 5, all plots contain the same event at E™ 1ss /-\/Ht ~ 
11 GeV 1//2 . The 'leptonic' backgrounds shown in the figures are those calculated from 
the Monte Carlo simulation, using the MC calculation of the cross section and normalized 
to 4.7 fb . The number of events observed in each of the six signal regions, as well as 
their Standard Model background expectations are shown in Table 3. Good agreement is 
observed between SM expectations and the data for all six signal regions. Table 3 also shows 
the 95% confidence-level upper bound -/Vg|j^ max on the number of events originating from 
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Signal region 


7j55 


8j55 


9j55 


6j80 


7j80 


8j80 


Multi-jets 


91±20 


10±3 


1.2±0.4 


67+12 


5.4+1.7 


0.42+0.16 


it -> q£, U 


55±18 


5.7±6.0 


0.70±0.72 


24+13 


2.8+1.8 


0.38+0.40 


W + jets 


18±11 


0.81±0.72 


0+0.13 


13+10 


0.34+0.21 


0+0.06 


Z + jets 


2.7±1.6 


0.05±0.19 


0+0.12 


2.7+2.9 


0.10+0.17 


0+0.13 


Total Standard Model 


167+34 


17±7 


1.9+0.8 


107+21 


8.6+2.5 


0.80+0.45 


Data 


154 


22 


3 


106 


15 


1 


KsM,max ( ex P) 
^B 5 SM, max (Obs) 


72 
64 


16 
20 


4.5 
5.7 


46 
46 


8.4 
15 


3.5 
3.8 


S, mB • A • e (exp) [ft] 
<?bsm^-A-£ (obs) [ft] 


15 
14 


3.4 
4.2 


0.96 
1.2 


9.8 
9.8 


1.8 
3.2 


0.74 
0.81 




0.64 


0.27 


0.28 


0.52 


0.07 


0.43 



Table 3. Results for each of the six signal regions for an integrated luminosity of 4.7 fb -1 . The 
expected numbers of Standard Model events are given for each of the following sources: multi-jet 
(including fully hadronic it), semi- and fully- leptonic it decays combined, and W and Z bosons 
(separately) in association with jets, as well as the total Standard Model expectation. The un- 
certainties on the predictions show the combination of the statistical and systematic components. 
Where small event counts in control regions have not made it possible to determine a central value 
for the expectation, an asymmetric bound is given instead. The numbers of observed events are also 
shown. The final five rows show the statistical quantities described in the text. Both the expected 
(exp) and the observed (obs) values are shown for max and cr||^j max x A x e. 

sources other than the Standard Model, the corresponding upper limit Cgg^ max X A X e on 
the cross section times efficiency within acceptance (which equals the limit on the observed 
number of signal events divided by the luminosity) and the p- value for the Standard-Model- 
only hypothesis (psm)- 

In the absence of significant discrepancies, limits are set in the context of two super- 
symmetric (SUSY) models. The first is the tan/3 = 10, Aq = and /i > slice of the 
MSUGRA/CMSSM parameter space. The second is a simplified SUSY model with only a 
gluino octet and a neutralino Xi within kinematic reach. Theoretical uncertainties on the 
SUSY signals are estimated as described in Section 5. Combined experimental systematic 
uncertainties on the signal yield from jet energy scale, resolution, and event cleaning are 
approximately 25%. Acceptance times efficiency values are tabulated elsewhere [49]. 

The limit for each signal region is obtained by comparing the observed event count with 
that expected from Standard Model background plus SUSY signal processes, taking into 
account all uncertainties on the Standard Model expectation, including those which are 
correlated between signal and background (for instance jet energy scale uncertainties) and 
all but theoretical cross section uncertainties (PDF and scale) on the signal expectation. 
The combined exclusion regions are obtained using the CL S prescription [50], taking the 
signal region with the best expected limit at each point in parameter space. The 95% confi- 
dence level (CL) exclusion in the tan (3 = 10, A = and fj, > slice of MSUGRA/CMSSM 
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is shown in Figure 7. The Her band surrounding the expected limit shows the variation 
anticipated from statistical fluctuations and systematic uncertainties on SM and signal 
processes. The uncertainties on the supersymmetric signal cross section from PDFs and 
higher-order terms are calculated as described in Section 5, and the resulting signal cross 
section uncertainty is represented by ±1<j lines on either side of the observed limit. 7 

The analysis substantially extends the previous exclusion limits [13, 16, 17] for mo > 
500 GeV. For large mo, the analysis becomes independent of the squark mass, and the 
lower bound on the gluino mass is extended to almost 840 GeV for large rrig. 8 In the 
simplified model gluinos are pair-produced and decay with unit probability to t + i+Xi. 
In this context, the 95% CL exclusion bound on the gluino mass is 870 GeV for neutralino 
masses up to 100 GeV. 

9 Summary 

A search for new physics is presented using final states containing large jet multiplicities 
in association with missing transverse momentum. The search uses the full 2011 pp LHC 
data set taken at y/s = 7 TeV, collected with the ATLAS detector, which corresponds to 
an integrated luminosity of 4.7 fb -1 . 

Six non-exclusive signal regions are defined. The first three require at least seven, eight 
or nine jets, with pr > 55 GeV; the latter three require at least six, seven or eight jets, 
with > 80 GeV. In all cases the events are required to satisfy E™ 1ss /\JHt > 4 GeV 1 / 2 , 
and to contain no isolated high-px electrons or muons. Investigations on the enlarged data 
sample have resulted in improvements compared to a previous measurement using a similar 
strategy. In particular, inclusion of events with smaller jet-jet separation increases the 
acceptance for signal models of interest by a factor two to five, without without significantly 
increasing the systematic uncertainty. 

The Standard Model multi-jet background is determined using a template-based method 
that exploits the invariance of E™ 1ss /\/Ht under changes in jet multiplicity, cross-checked 
with a jet-smearing method that uses well reconstructed multi-jet seed events from data. 
The other significant backgrounds — tt + jets, W + jets and Z + jets — are determined 
using a combination of data-driven and Monte Carlo-based methods. 

In each of the six signal regions, agreement is found between the Standard Model pre- 
diction and the data. In the absence of significant discrepancies, the results are interpreted 
as limits in the context of i?-parity conserving supersymmetry. Exclusion limits are shown 
for MSUGRA/CMSSM, for which, for large mo, gluino masses smaller than 840 GeV are 
excluded at the 95% confidence level. For a simplified supersymmetric model in which both 
of the pair-produced gluinos decay via the process g — > t + t + Xi, gluino masses smaller 
than about 870 GeV are similarly excluded for Xi masses up to 100 GeV. 



7 Previous analyses have a slightly different presentation of the effect of the signal cross section uncer- 
tainty. In Refs. [13, 16, 17] the effect of the signal cross section uncertainty was folded into the displayed 
limits and so was not shown separately. 

s Limits on sparticle masses quoted in the text are those from the lower edge of the 1 a signal cross 
section band rather than the central value of the observed limit, so can be considered conservative. 
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Figure 2. Jet multiplicity distributions for the tt + jets validation regions (left) and control 
regions (right) before any jet multiplicity requirements, for a jet pt threshold of 45 GeV (top), 
55 GeV (middle) and 80 GeV (bottom). 



- 14 - 



I L dt ~ 4.7 fb 
ATLAS 



Data 2011 (\!s= 7TeV) 
Total SM prediction 
Alpgen W-> ftv 1 
Alpgen tf-> qljl 
Alpgen Z-> fee,fi|i,ti) 1 
Alpgen W-> (e,i)v 
SUSY m 0= 2960, m„ 2= 240E 




W validation region 
p T > 55 GeV jets 











— • - • • 1 









1 23456789 10 

Number of jets 



1 10 7 

iSio 6 

10 5 
10 4 
10 3 
10 2 
10 

1 

10 1 



J L dt ~ 4.7 fb 
ATLAS 



I i i i i I i i i i I i i i i I i i i i I i 
Data 2011 (\ls=7TeV) ' 
— Total SM prediction 
| Alpgen W-> ftv 

□ Alpgen tt-> qljl 

| Alpgen Z—> (ee,fi|i,rr) 

^ Alpgen W^> (e,t)v 

.... SUSY m 0= 2960, m ])2= 240| 

W control region 
p T > 55 GeV jets 



7 8 9 10 
Number of jets 



(a) 



(b) 



J L dt ~ 4.7 fb 
ATLAS 




Data 2011 (\ls = 7TeV) j 
— Total SM prediction 

| Alpgen W-> 
■ Alpgen tt^ qljl I 
| Alpgen Z-> (ee,n&l,TC) 
B Alpgen W-> (e,i)v 1 
.... SUSY m 0= 2960, m 1B =240: 

W validation region 
p T > 80 GeV jets [ 



7 8 9 10 
Number of jets 




Data 2011 (\Js=7TeV) 
Total SM prediction 
Alpgen ftv 

Alpgen tt-> ql,ll 
Alpgen Z-> (ee,^fi,Tt) 
Alpgen W-> (e,i)v 
SUSY m 0= 2960, m„ 2 =240: 

W control region 

80 GeV jets 



7 8 9 10 
Number of jets 



(c) 



(d) 



Figure 3. Jet multiplicity distributions for the + jets validation regions (left) and control 
regions (right) before any jet multiplicity requirements, and for a jet pr threshold of 55 GeV (top) 
and 80 GeV (bottom). 
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Figure 4. As for Figure 3 but for the Z + jets validation regions and control regions. 
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Figure 5. The distribution of the variable E™ lss j \J Ht for each of the six different signal regions 
defined in Table 1, prior to the final Ef iss /^/Hr > 4 GeV 1/2 requirement. 
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Figure 6. The distribution of jet multiplicity for jets with > 55 GeV (a) and those with 
p T > 80 GeV (b). Only events with E T liss /y/Hr > 4 GeV 1/2 are shown. 
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Figure 7. Combined 95% CL exclusion curves for the tan j3 = 10, Aq = and /.i > slice of 
MSUGRA/CMSSM (a) and for the simplified gluino-neutralino model (b). The dashed grey and 
solid red lines show the 95% CL expected and observed limits respectively, including all uncertainties 
except the theoretical signal cross section uncertainty (PDF and scale). The shaded yellow band 
around the expected limit shows its ±1<t range. The ±lcr lines around the observed limit represent 
the result produced when moving the signal cross section by ±ler (as defined by the PDF and 
scale uncertainties). The contours on the MSUGRA/CMSSM model show values of the mass of 
the gluino and the mean mass of the squarks in the first two generations. Exclusion limits are also 
shown from previous ATLAS searches with >2, 3 or 4 jets plus _E™ 1SS [16], multi-jets plus E™ lss [13] 
or with same-sign dileptons [46] and from LEP [47] in (a) . The lower plot shows limits from ATLAS 
searches with same-sign dileptons [46] or with one-lepton plus fe-jet [48]. 
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A Event displays 

A display of an event that passes the 9j55 and 7j80 signal region selections can be found 
in Figure 8. A display of an event that passes all signal region selections can be found in 
Figure 9. 
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(a) 



Figure 8. A display of an event which passes the 9j55 and 7j80 signal region selections. The 
event has E^ iss /^/E^ of 4.1 GeV 1/2 , H T of 1.47 TeV and E^ iss of 157 GeV. 
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(a) 



Figure 9. A display of an event which passes all signal region selections. The event has E!f lss / \/Ht 
of 11.6 GeV 1/2 , H T of 1.17 TeV and E^ iss of 397 GeV. One of the jets, with p T of 107 GeV is b 
tagged. The event also contains a muon with px of 90 GeV, overlapping with a jet. 
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